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ABSTRACT 


Turbulent  characteristics  of  the  first 
62  meters  of  the  atmosphere  in  the  vicinity  of 
the  U.  S.  Army  Electronics  Research  and  Develop¬ 
ment  Activity's  Meteorological  Research  Tower 
are  established  for  neutral  conditions.  The 
assunption  was  made  that  the  roughness  length  is 
a  constant,  but  dependent  upon  wind  direction, 
fetch,  and  the  he'C'.c  of  the  roughness  elements. 

Data  are  presented  for  five  recording 
periods  during  the  late  winter  and  early  spring 
of  1958  and  1961,  Computations  of  the  basic  wind 
profile  and  turbulence  parameters  are  presented 
in  tabular  form, 
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INTRODUCTION 


The  notion  of  the  atmosphere  near  the  surface  of  the  earth  is  rarely  as 
simple  as  described  by  the  theory  of  laminar  flow.  Normally,  flow  in  the 
atmosphere  consists  of  a  fairly  simple  mean  motion  on  which  is  superimposed 
extremely  complicated  secondary  or  eddy  motions  of  oscillatory,  but  not  nec¬ 
essarily  periodic,  character.  The  superimposed  eddy  motion  results  in  the 
often  chaotic  wind  profile  and  turbulent  spectrum  observed  in  the  boundary 
layer. 

Fluctuations  about  the  mean  flow  in  the  atmosphere  are  the  function  of 
a  number  of  conditions,  including  the  air-earth  interface  and  the  exchange 
of  energy  between  the  atmosphere  and  the  underlying  surface.  These  boundary 
conditions  plus  the  prevailing  lapse  rate  are  the  determining  factors  for 
variations  in  the  local  wind  profile.  The  basis  for  any  study  concerning 
turbulent,  or  eddy  motion  in  the  atmosphere  is  the  establishment  of  initial 
boundary  conditions.  This  report  is  concerned  with  establishing  these  basic 
parameters  for  the  immediate  vicinity  of  the  U.  S.  Army  Electronics  Research 
and  Development  Activity's  (USA  EROA)  Meteorological  Research  Tower. 


TOE  METEOROLOGICAL  TOWER  FACILITY 


The  meteorological  tower  facility  was  established  for  basic  research 
with  respect  to  wind  effects  on  unguided  rockets,  including  wind  profile 
studies,  the  spectra#  of  turbulence,  and  the  prediction  of  the  local  wind 
profile.  The  tower  instrumentation,  the  surrounding  torvain,  and  the 
climatology  of  White  Sands  Missile  Range  (WSMR) ,  Now  Mexico,  have  been 
described  adequately  elsewhere  [1,  2], 


WIND  PROFILE  HYPOTHESES 


The  two  hypotheses  most  conmonly  used  to  describe  the  wind  profile  in 
the  boundary  layer  are  the  power  law  and  the  logarithmic  law.  The  pewer 
law,  better  known  as  the  "seventh  root  law"  or  Schmidt's  conjugate  power 
law,  normally  takes  the  form 


(-1 r)P 

where  u  is  the  mean  wind  speed  at  height  z, 


Uj  is  the  mean  wind  speed  at  height  zj,  and 
p  is  the  profile  index. 
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Its  use  is  limited  to  profiles  over  surfaces  which  are  either  bare  or 
covered  with  very  short  vegetation,  The  profile  index  p  is  also  dependent 
upon  the  vertical  temperature  gradient, 

The  logarithmic  profile  is  based  upon  a  steady  two-dimensional  motion 
in  the  proximity  of  a  surface  z  »  0  with  H  as  a  function  of  z  only.  In  fully 
developed  turbulent  flow,  tho  object  is  to  find  a  velocity  profile  u  (z) 
which  is  consistent  with  a  given  shear  stress  and  under  the  assumption  that 
the  temperature  gradient  is  near  adiabatic.  The  equation  has  taken  many  forms, 
the  most  common  being  attrihuted  to  Rossby  [3]: 


for  a  fully  rough  surface,  where 

IT  ■  mean  wind  speed, 

* 

»  friction  velocity, 


k  «  von  Karman ’s  constant, 
z  *  height, 

z0  »  roughness  length. 

Utilizing  the  power  and  logarithmic  laws  in  the  preceding  forms,  ex¬ 
pressions  for  the  boundary  layer  parameters  can  be  derived. 

The  Roughness  Length:  The  roughness  length,  z_,  a  constant  of  integration, 
is  related  to  the  neigh t  of  the  surface  irregularities  by  z0  :  e/30,  This 
implies  that  the  mean  velocity  tends  toward  zero  at  a  height  depending  upon 
tho  average  length  of  the  roughness  elements.  The  roughness  length  z0  can 
be  easily  determined  from  the  IT,  log  z  plot  under  adiabatic  or  neutral  lapse 
conditions. 


The  Surface  Shearing  Stress ,  Friction  Velocity,  and  Prag  Coefficient;  The 
keynoi'ds  stresses  indicate  that  fluctuations  in'  velocity  cause  transport  of 
momentum  arms*  a  surface  in  a  fluid,  Ir  general,  the  Reynolds  stresses  out¬ 
weigh  the  viscous  stresses,  which  often  may  be  neglected  in  problems  of  turbu¬ 
lent  motion.  The  horizontal  shearing  stress  [3]  is  given  by: 

i  ■  -pu'w' 


where  p  is  mean  density  and  u’w'  are  mean  fluctuating  or  eddy  velocities, 

The  shearing  stress  is  assumed  to  be  independent  of  height  within  the  first 
50-100  meters  above  tho  surface,  according  to  Lettau  [4],  and  related  to  the 
friction  velocity  by 


(luTH)* 
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and  to  the  surface  shearing  stress  and  the  coefficient  of  drag  by 


T  “  T0  +  2 

ax 

T0  =  £  Pqfuf 

where  ug  is  the  wind  velocity  "near"  the  surface.  Since 

CD  =  2jo  , 

pH* 

it  follows  that 


assuming  the  Uossby  profile. 

The  Eddy  Viscosity;  Sutton  [5]  states  that  a  basic  step  in  turbulence  research 
is  to  adopt  the'  fundamental  ideas  of  the  kinetic  theory  of  gases  by  expressing 
the  transfer  of  momentum  or  any  other  suitable  entity,  by  means  of  virtual 
coefficients  of  viscosity,  conductivity,  and  diffusivity,  defined  in  much  the 
spne  way  as  their  molecular  counterparts,  These  exchange  coefficients,  ex¬ 
pressions  for  the  turbulent  flux,  may  be  derived  independently  of  any  theory 
of  the  structure  of  eddy  motion.  The  eddy  viscosity,  or  the  turbulent  trans¬ 
fer  of  momentum  by  eddies  giving  rise  to  an  interna]  fluid  friction,  is  defined 
as 


I  3  O  +  Km)  dO  -  KM  du  if  v  «Km 
P  dz  uz 


du 

~rr 


where  Km  is  the  eddy  viscosity.  From  the  above  relationship  and  the  .tosshy 
profile,  the  edd.-  v.'  siosity  car  br  drfin:d  .-.s 

kM  “  k"*' 

which  leads  to  a  valid  approximation  for  boundary  conditions. 


DATA  ANALYSIS 


Wind  and  temperature  profiles  originally  observed  for  spectral  analysis 
purposes  at  the  tower  site  were  used  to  determine  boundary  conditions  and  to 
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investigate  the  application  of  the  p  ro :  1  i  c  hypotheses  to  the  haste  data 
The  raw  data  we  it  recorded  on  five  different  days,  and  the  length  of  re¬ 
cording  varied  frc..n  or.e  to  six  hours  One-hour  samples  were  obtained  on 
12  February  19SS,  13  Februaiy  1958,  and  7  March  1958,  and  consisted  of 
wind  velocity  and  temperatures  for  nine  tower  levels  Six-hour  samples 
were  obtained  on  IS  and  16  March  19b!  for  four  tower  levels.  The  1958 
data  were  reduced  as  five-second  visual  averages  for  the  wind  velocities, 
with  temperatures  read  every  44  seconds.  The  1961  data  were  .educed  as 
four-second  and  sixty-second  averages  for  the  winds  and  temperatures, 
respectively, 

All  data  were  punched  on  Hollerith  cards  for  machine  analysis-  Wind 
and  Temperature  profiles  by  the  hour  and  by  ten-minute  intervals  were  then 
computed  using  the  Philco  200(1  computer  The  fifteen  hourly  means  and  de¬ 

viations  are  presented  in  Tables  I  to  V.  The  90  ten-minute  mean  profiles 
were  used  to  establish  boundary  conditions  for  the  tower  data  (Tables  VI 
to  XI  , 


BOUNDARY  CONDITIONS 


To  establish  the  basic  turbulence  parameters  for  the  tower  area,  the 
wind  and  temperature  profiles  were  inspected  for  adiabatic  or  neutral  con¬ 
ditions  Of  the  90  profiles  available,  11  were  found  thnt  mot  the  criteria, 
The  11  wind  profiles  fere  plotted  (figures  1-11)  on  semi -logarithmic  paper, 
and  the  line  of  regression  determined  by  a  least  squares  fit.  The  IT,  log  z 
intercept  was  taken  as  the  roughness  length  z0  The  equivalent  surface 
roughness  vs  wind  speed  at  4  3  and  4,6  meters  is  shown  in  Figure  12, 

The  average  height  o 7  the  roughness  elements  surrounding  the  tower 
facility  is  approximately  three  meters,  giving  a  theoretical  z0  of  ten 
centimeters,  which  is  in  close  agreement  with  the  observed  data  and  with 
the  results  of  other  experimenters  in  the  field 

Assuming  the  roughness  length  to  be  valid  for  the  tower  area,  values 
of  z0  picked  from  rhe  curve  (Figure  12)  were  used  to  compute  the  Friction 
velocities  for  each  ten -minute  profile  The  remainder  of  the  boundary 
parameters  under  discussion  were  also  computed  at  this  time  and  are  pre¬ 
sented  in  Tables  XI  through  XV.  All  calculations  were  for  a  height  range 
of  4  3  to  4  6  meters  above  the  surface  and  a  fetch’  of  40  to  70  meters 


DISCUSSION 


A  careful  examination  of  the  boundary  parameter.-,  indicates  a  dependence 
on  the  mean  wind  speed;  however,  in  the  case  of  the  coefricient  of  drag,  Cp, 
it  is  apparent  that  the  wind  direction  has  a  decided  effect.  The  larger 
values  of  Cjj  0,029  to  0.032)  tend  to  occur  when  the  mean  wind  direction  is 
from  the  quadrant  in  which  the  largest  roughness  elements  are  with  respect 
to  the  tower  (Figure  13i 


’Fetch  is  defined  as  the  upwind  distance  from  nearest  obstruction  to  the 
ins  t  rumen  tat  i  on 
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Remarks :  Wind  speeds  in  m  sec  ;  Temperatures  in  degrees 
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Reitarks:  Wind  speed  in  m  sec-1;  Temperature  in  degrees 
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Remarks :  Wind  speeds  in  m  sec”  •  Temperatures  in  degrees 
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Remarks:  Wind  speeds  in  a  sec  .  Temperatures  in  degrees 
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1221-1230  874.5  1029.1  6.3  45  0.018  2.05  0.828  283.5 
1231-1240  874.2  1028.0  5-3  3?  0.017  1.40  0.681  196.1 
1241-1250  873.9  1027.6  6.6  48  0.019  2.34  O.885  316.8 
1251-1300  873.6  1027.3  6.0  43  0.018  1.87  0.791  258.0 
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2321-2330  871.8  1081.4  7-5  54  0.020  3.15  0.929  405.0 
2331-2340  871.8  108l.4  8.2  58  0.021  3.69  0.998  475.6 
2341-2330  871.8  1080.6  9.5  66  0.022  4.72  1.135  627-0 
2351-2400  871.8  1080.6  9.0  64  0.022  4.42  1.101  576.0 
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Figure 
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WIND  SPEED  m/'sec 

LOGARITHMIC  WIND  PROFILE.  1364  -  1403  MST 
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10GARI JHH  Sf  HEIGHT.  Meters 


V  I H  0  SPEED  m / S  8  c 

LOGARITHMIC  WIND  PROFILE  FOR  1404  -  1413  MST 
7  MARCH  1058,  y  -  0.  I498b<>*,‘207,‘ 

Figure  3 
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WIND  SPEED  n/set 

LOGARITHMIC  WIND  PROFILE  FDR  1424  -  1433  MST 
7  MARCH  1058.  y  -  0.2484e°'*861“ 

Figure  4 
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Figure  9 
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250 

DIRECTION  IN  DEGREES 
V I ND  DIRECTION  vs  COEFFICIENT  OF  DRAG 
Figure  13 


All  values  of  roughness  length,  friction  velocity,  drag  coefficient, 
surface  stress  and  eddy  viscosity  presented  in  Tables  XI  to  XV  appear  to 
be  reasonable  and  are  in  good  agreement  with  those  of  other  investigators, 
notably  Deacon  [6]  and  Sheppard  [7]  who  found  roughness  lengths  of  9  cm 
for  thick  grass  up  to  SO  cm  high,  while  drag  coefficients  for  thick  grass, 
assuming  the  Rossby  {3|  profile,  were  found  to  be  on  the  order  of  0.032, 

The  Wind  Profile:  Since  neither  the  logarithmic  profile  nor  the  power  law 
is  actually  valid  under  diabatic  conditions,  few  conclusions  may  be  drawn 
concerning  the  observed  wind  profiles,  It  can  be  generally  stated  that  the 
power  law  would  provide  the  best  fit  for  this  particular  sample  of  data, 
which  is  in  agreement  with  Swanson  and  lloidale  [8], 

Test  for  a  Fully •  Rough  Surface;  With  the  concept  of  fully- rough  flow,  the 
influence'  o(  viscosity  fs  'considered  to  be  negligible,  Nikuradse  [9]  has 
proposed  tests  for  smooth  and  rough  surfaces,  which  can  be  stated  as 


Smooth  Flow:  J^»2o  <  0,13 
v 


Fully- rough  Flow;  '‘*2o  >2.5 
v 

assuming  that  the  relation  z0  «  c/30  can  be  accepted. 

The  quantity  u*zQ  is  known  as  the  macro- viscosity  and  is  denoted  by  N, 
Nikuradse’s  criteria  can  thus  be  stated 


Smooth  Flow:  N  <  0.13  v  ;  0.02  cm2  see-' 
Rough  Flow:  N  >  2,5  v  '  0,4  cm2  sec-1. 


Values  of  N  are  presented  in  Tables  XI  through  XV. 

It  is  seen  that  all  data  meet  the  rough  flow  criteria  and  the  hypothesis 
assumed  is  valid 


Turbulent  Flow  Indicators.  Cramer  [10]  shows  that  the  standard  deviation  of 
wind  "liliection  is  a  gootf  indicator  of  stability  and,  of  course,  turbulent 
intensity.  Cramer’s  estimates  are  presented  in  Table  XVI, 
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TABLE  XVI 


ESTIMATED  RANGE  IN  STANDARD  DEVIATION  OF  AZIMUTH  KIND  DIRECTION,  °A, 
NEAR  GROUND  LEVEL  FOR  VARIOUS  STABILITY  STRATIFICATIONS 


St  jo  itity 
Stratification 


Smooth  Site 


Rough  Site 


Extremely  stable 

2- 

4 

Moderately  stable 

4- 

6 

Near  Neutral 

6- 

8 

Moderately  unstable 

10- 

15 

Extremely  unstable 

20- 

25 

2-  6 
7-15 
10- IS 
15*20 
25-30 


It  is  clear  from  Tables  I  to  X  and  Figure  14  that  the  data  agree  well 
with  Cramer's  hypothesis,  This  is  reflected  by  the  gustiness  ratio  or 
intensity  of  turbulence 


I 


°TT 


u 


where  I  is  the  gustiness  ratio  and  8u  is  the  standard  deviation  about  the 
mean  wind  speed  tt,  Figure  14  is  a  plot  of  the  standard  deviation  of  direc¬ 
tion  vs  stall!  lity,  The  stability  ratio  represents  the  ratio  of  tempsrature 
difference  "n  degrees  Celcius  between  4.3  or  4.6  meters  and  25,6  or  31,7 
meters  to  the  square  of  the  mean  wind  speed  at  19,3  or  19.5  meters,  depending 
upon  the  data  sample  used. 

Another  indicator  of  turbulent  motion  in  the  boundary  layer  is  an  ap¬ 
parent  inverse  relationship  between  the  variance  of  wind  direction  and  the 
eddy  visensitv  Tt  appears  that  as  the  variance  increases,  the  eddy  vis¬ 
cosity  decreases,  indicating  a  breakdown  of  the  viscous  stresses  with  in¬ 
creased  mi xing. 

.^rtVhile  boundary  conditions,  arc  dependent  upon  wind  speed  and  are  neces¬ 
sary  for  establishing  the  characteristics  of  turbulence  for  an  area,  it  is 
obvious  that  the  variance  of  wind  direction  is  the  determining  factor  in 
turbulent  intensities  and  the  propagation  of  turbulent  motion  downstream, 
Bhckadar,  Panofsky,  et  al.  [II]  show  from  Taylor's  Hypothesis  that 
Lagrangian  and  Euleriun  standard  deviations  of  azimuth  are  approximately 
equal  during  daytime,  and  that  measurements  at  a  fixed  site  are  valid 
some  distance  downstream,  Thus,  Taylor's  Hypothesis  can  be  used  to  evaluate 
the  behavior  of  the  wind  after  it  has  passed  a  fixed  point. 
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STA3ILITY  RATIO  T  -  T  VS.  STANDARD  DEVIATION  OF  Dl SECTION 


CONCLUSIONS 


TUrbulent  characteristics  of  the  atmospheric  boundary  layer  in  the 
vicinity  of  the  U,  S.  Army  Electronics  Research  and  Development  Activity 
Meteorological  Research  Tower  are  comparable  to  those  observed  at  other 
locations.  It  can  be  concluded  that  the  roughness  length  increases  with 
increasing  wind  speed  and  varies  with  wind  direction.  This  is  attributed 
to  changes  in  the  fetch  and  the  height  of  the  roughness  elements  from  dif¬ 
ferent  exposures  due  to  wind  direction  shifts.  The  coefficient  of  drag  of 
the  surface  varies  with  wind  direction,  exhibiting  the  same  characteristics 
as  the  roughness  length  with  relation  to  the  fetch  and  roughness  elements. 

The  standard  deviation  of  wind  direction  was  found  to  vary  with  the 
lapse  rate  while  the  eddy  viscosity  appears  to  be  inversely  proportional 
to  all  standard  deviation  of  direction.  All  boundary  -.onditions  were 
found  to  be  dependent  upon  wind  speed  and  to  some  degree  dependent  on 
wind  direction.  The  standard  deviation  of  wind  direction  was  determined 
to  be  a  good  indicator  of  tuibulence. 
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LIST  OF  SYMBOLS 


A  Wind  Direction,  Degrees 

Me  on  Wind  Direction,  Degrees 
Cjj  Coefficient  of  Drag 

I  Gustinoss  Ratio 

Eddy  Viscosity 
N  Macro  Viscosity 

P  Ambient  Pressure 

T  Temperature,  Degrees  Celsius 

k  Karman’s  Constant 

p  Wind  Profile  Index 

u  Horizontal  Wind  Speed 

u*  Instantaneous  Horizontal  Velocity 

u  Mean  Wind  Speed 

u  Friction  Velocity 

w  Vertical  Velocity 

w’  Instantaneous  Vertical  Velocity 

z  Height 

z0  Roughness  Length 

e  Roughness  Element 

p  Density 

t  Shearing  Stress 

t0  Surface  Shearing  Stress 

v  Kinematic  Viscosity 
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